Loss of molecules in magneto-electrostatic traps due to nonadiabatic transitions 
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We analyze the dynamics of a paramagnetic, dipolar molecule in a generic "magneto-electrostatic" 
trap where both magnetic and electric fields may be present. The potential energy that governs 
the dynamics of the molecules is found using a reduced molecular model that incorporates the main 
features of the system. We discuss the shape of the trapping potentials for different field geometries, 
as well as the possibility of nonadiabatic transitions to untrapped states, i.e., the analog of Majorana 
transitions in a quadrupole magnetic atomic trap. Maximizing the lifetime of molecules in a trap 
is of great concern in current experiments, and we assess the effect of nonadiabatic transitions on 
obtainable trap lifetimes. 

PACS numbers: 



I. INTRODUCTION 

A primary tool for the study of ultracold matter is the 
magnetic trap, which can confine paramagnetic atoms on 
time scales of seconds. Using the techniques of laser cool- 
ing, atoms can be loaded into magnetic traps — such as 
the quadrupole trap created with anti-Helmholtz coils — 
at cold enough temperatures so that only modest mag- 
netic field gradients are required to confine the atoms. 
Once loaded, the lifetime in a quadrupole trap is typi- 
cally limited by collisions, thermal background radiation, 
or Majorana transitions [1], wherein atomic spins pass- 
ing near the zero in the field minimum are not able to 
adiabatically follow the field direction. To such an atom, 
the quantization axis is lost and it may emerge from the 
trap zero "spin-flipped" to an untrapped state [2]. 

The Majorana spin-fiip loss mechanism can be severe 
enough to quench trap lifetime before forced evapora- 
tive cooling can bring the atomic system to quantum de- 
generacy. Such cooling is thought to be a necessary re- 
quirement to reach the strongly correlated regimes that 
hold much promise for investigating novel ultracold col- 
lisions, chemistry, connections to condensed matter, and 
quantum information processing with particles beyond 
presently studied atoms, e.g., polar molecules [3]. To 
prevent such loss in atomic systems, evaporative cool- 
ing is instead performed in a secondary magnetic trap 
of the loffe-Pritchard or "baseball" [2, 4] configuration 
which does not possess a vanishing field magnitude near 
the trap minimum, and thereby preserves a quantization 
axis for all particles confined in the trap. Additional tech- 
niques to minimize Majorana loss include moving trap 
minima in time so that the atoms experience a non-zero 
field on average [5] , and confining in storage rings wherein 
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the centrifugal barrier shifts the trap center away from 
field minima [6]). 

Large samples of molecules in their absolute rovibronic 
ground state have yet to be produced at temperatures 
below 10 mK [3]. Sample confinement therefore re- 
quires very large magnetic fields, even for the case of 
a quadrupole trap [7]. Electrostatic [8, 9] and mag- 
netostatic [10, 11] traps have successfully been demon- 
strated for molecules, and inhomogeneous static electric 
and magnetic fields have been combined to increase trap 
depth [7] . Lifetimes have so far been limited by collisions 
or blackbody radiation, but once evaporative or other 
cooling schemes [3, 12] are employed, Majorana tran- 
sitions will be a dominant cause of loss in all the cur- 
rently employed traps for ground state molecules. Since 
the large fields necessary to create, e.g., loffe-Pritchard 
traps, are generally unobtainable, we must carefully in- 
vestigate the lifetime limits imposed by Majorana tran- 
sitions in the currently realizable quadrupole traps and 
explore possible schemes to avoid such loss mechanisms 
for real molecular systems. 

For cold molecules that possess both magnetic and 
electric dipole moments, novel traps formed from both in- 
homogeneous magnetic and electric fields may be utilized 
to confine molecules in ways that mitigate the unwanted 
Majorana transitions. In this paper we explore the ba- 
sic physics of traps for molecules in which both magnetic 
and electric fields are present, which we show can mit- 
igate Majorana loss even for simple field configurations 
consisting of superimposed quadrupole and homogenous 
fields. Using a reduced, but analytically solvable, model 
of a molecule in both electric and magnetic fields, we ex- 
plore various trap geometries that can be generated from 
common coil and electrode configurations. We identify 
the zone in the trapping potential where nonadiabatic, 
Majorana-like transitions can occur and we introduce ap- 
proximate expressions to assess such loss. More broadly, 
we note that these results are of interest to the design 
of Stark decelerators [13], where slowed molecules expe- 
rience rapidly changing field configurations. Our main 
focus will be on diatomic molecules governed by Hund's 
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coupling case (a) . since in this case the electric and mag- 
netic dipole moments are nontrivially related. In particu- 
lar, we use the OH molecule as a representative molecular 
system, but we briefly address case (b) molecules as well. 

II. MOLECULAR STRUCTURE IN THE 
PRESENCE OF CROSSED FIELDS 

Before discussing the Majorana spin-flip problem for 
representative traps and molecules, we first examine per- 
turbations to the molecular structure in crossed electric 
and magnetic fields. We consider a diatomic, heteronu- 
clear molecule that has a permanent electric dipole mo- 
ment fip lying along its molecular axis. The molecule 
is moreover assumed to be paramagnetic, with magnetic 
dipole moment /Xj„. In the presence of spatially vary- 
ing magnetic (jB) and electric {£) fields, the molecule is 
governed by the Hamiltonian 

H^Ho--fim-B{r)^fi,-eir), (1) 

where Hq describes the internal workings of the molecule, 
including lambda doubling, fine structure, and hyperfine 
structure, if applicable. The eigenencrgies of H will vary 
in space according the spatial variation of the fields, and 
these varying energies define the trap potentials. 

Eigenvalues of H can be constructed to any desired 
degree of accuracy, generating realistic trap potentials 
for any desired molecule. Such potentials were generated 
for OH in Ref. [7] according to this procedure. Presently, 
however, we are interested in the general features of the 
traps, and so will opt for a simplified molecular structure, 
albeit one which encapsulates all the relevant physics. We 
will explicitly consider Hund's case (a) molecules with ^11 
electronic symmetry (as in OH) as well as Hund's case (b) 
molecules with symmetry. These electronic structures 
are the most commonly encountered examples of open- 
shell diatomics, and throughout we will neglect hyperfine 
structure. 

A. Hund's case (a) ^11 molecule 

Case (a) molecules are characterized by a strong spin- 
orbit interaction that couples the electronic spin to the 
electronic orbital angular momentum, and thus also to 
the molecular axis. The relevant quantum numbers for a 
case (a) molecule are \ {AT,)JMfl), where A and S are the 



projections of the orbital and spin angular momenta of 
the electrons on the molecular axis, respectively; J is the 
total angular momentum of the molecule, resulting from 
the addition of the orbital and spin angular momentum 
of the electrons, and the rotation of the nuclei; M is the 
projection of J on the laboratory axis; and O its projec- 
tion on the molecular axis. We will suppress A and E 
quantum numbers which take assumed values hereafter. 

In the absence of an external electric field, the molecule 
is also an eigenstate of parity, characterized by another 
quantum number, e = ±1, which is conventionally de- 
noted by e and /, respectively: 

\{e/f)JMn) = l/V2{\JMn) + e\JM - n)), 

where Q = \Q\ is defined to be positive and the parity 
of the state is given by e(— 1)'^^*. The energy gap be- 
tween these two states due to A-doubling [14] is denoted 
by A (which equals 27r x 1.7 GHz for OH). OH molecules 
in their ground state are easily polarizable due to their 
relatively small A-doubling and large electric dipole mo- 
ment (/Ze = 1.67 D). In an electric field, the superposition 
of opposite parity states gives rise to an effective polar- 
ization. 

For simplicity, we will assume that J = 1/2 and thus 
fl = 1/2 (we will suppress them also in the notation), 
and will consider only the four states \f,M = -1-1/2), 
|/,M = -1/2), |e,M = +1/2) and \e,M = -1/2). We 
will further assert that the fields applied are small enough 
that J remains approximately conserved. Unfortunately, 
there exists an almost complete cancelation of orbital and 
spin magnetic dipole moment for a ^Ili/2 molecule, lead- 
ing to a nearly vanishing Zeeman interation. This defi- 
ciency in our model could be easily eliminated by con- 
sidering a ^n3/2 state, but at the expense of introducing 
additional M levels. To maintain analytical simplicity, 
we will formally use a ^lli/2 molecule, but ascribe to 
it a molecular g- factor of 4/5 which is the value corre- 
sponding to OH's ground state. This ad hoc parameter 
adjustment will not change the qualitative conclusions of 
the model, as we have checked by comparison with the 
results for a more complicated ^113/2 model. 

Using this molecular basis and without losing general- 
ity, we begin by assuming there is a magnetic field along 
the Z-axis at all points in space and an electric field in 
the XZ-plane which makes an angle a with respect to 
the magnetic field. The matrix elements of H are easily 
calculated (see Appendix A): 



H 



f-A/2 + Um -Uecosa -C/eSina \ 

— A/2 — Um — f/gsina UeCosa 

— [/,. cos a —Ue sin a A/2 + U„i 

\ — i/gSina Ugcosa A/2 — UmJ 



(2) 



3 



where 

Ue = .9eMe|^l 

= Sm/^sl^l, (3) 

are characteristic electric and magnetic energies and [Ib 
is the Bohr magneton. Geometric factors and any depen- 
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FIG. 1: Stark (top panel) and Zeeman (bottom panel) en- 
ergies as a function of the respective field magnitude for the 
four-state model defined in the text, with parameters corre- 
sponding to a OH molecule (except for a ~1 /ib mag- 
netic moment). The states' energies are ±A/2 in the absence 
of fields. In the case of Stark effect, both the weak and strong 
field-seeking energies are doubly-degenerate. The Zeeman ef- 
fect breaks the degeneracy except when U,n = A/2. The 
matrix in Eq. (2) is diagonal and each line in the figure cor- 
responds to one |(e//),A'/ = ±1/2) basis vector. 



where we have assumed the existence of four different 
nondegenerate eigenvalues, and labeled them in increas- 
ing order of energy (i.e., "Ui" labels the lowest energy 
eigenvalue) . The associated eigenvectors will be denoted 
|1), |2) |3) and |4). In the limit where one field vanishes, 
this expression reduces to the familiar Stark and Zeeman 
effects (see Fig. 1). 

If the electric and magnetic fields are everywhere par- 
allel, then the total energy is simply the sum of the inde- 
pendent Stark and Zeeman energies. In general, however, 
we must account for the angle a between these fields. The 
eigenenergies are plotted versus a in Figure 2 for the four 
states. We assume Ue and Um values that correspond to 
our model OH molecule in fields of |5| = 5 x 10'^ V/cm 



dence on the particular quantum numbers of the state 
have been gathered into a global g-factor which can in- 
terpreted as the effective coupling with the field. Explicit 
expressions are given in Appendix A. This matrix can be 
diagonalized analytically at every point in space, yielding 
the following energy eigenvalues: 



(4) 
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FIG. 2: (Color online) Dependence of the eigenenergies of 
the four-state model on the angle a between the electric and 
magnetic field directions for fixed values of their magnitudes 
(see righthand drawing). 



and |B| = 5 X lO'' G. This figure illustrates that the states 
experience avoided crossings in the presence of non-zero 
A-doubling for a = n/2. As we will see in further ex- 
amples, such competition between electric and magnetic 
effects enrich the trap potentials. A semi-classical treat- 
ment of this situation is given in Appendix B. 

As an aside, we note interesting features of the energies 
in crossed fields which may be useful in contexts beyond 
trapping. Applying a constant magnetic field of magni- 
tude such that Um — A/2 (see Fig. 1) causes two oppo- 
site parity states to be degenerate [15], as can be seen by 
taking the zero-electric-field limit Ue —>■ in Eq. (4)). Be- 
cause of this degeneracy, the Stark effect is immediately 
linear in a nonzero electric field, regardless of the angle 
between fields. The slope of this linear dependence, how- 
ever, may depend on a. This may be useful for increasing 
the efficiency of Stark decelerators [16]. 

Conversely, a Zeeman effect that is quadratic in \B\ 
can be obtained by applying a constant perpendicular 
electric field which induces an avoided crossing between 
two states. Rearranging terms in Eq. (4) and substituting 
a = 7r/2 into the eigenvalues, one can obtain U3^2 = 
±^JUl {Um T A/2)2. The constant electric field has 
induced a quadratic Zeeman effect in the vicinity of Um = 
A/2. This effect could be used to guide molecules moving 
in two dimensions: in a region where a linearly varying 
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f/l,2,3,4 = ± V^(A/2)2 + Ul +Ul± 2Um [(A/2)2 + {Ue COS a)^]' 
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magnetic field exists, a constant perpendicular electric 
field, for instance, induces minima at points where the 
magnetic field has the value \B\ = A/2fj.Bgm (^1-5 kG 
for our model OH). 



B. Hund's case (b) molecules 

Molecules like CaH, the first molecular system trapped 

by the buffer gas cooling method [10], or YbF, whose 
trapping could be useful to measure the electric dipole 
moment of the electron [17], belong to the category of 
molecules known as Hund's case (b) [14]. We consider 
here a pure case (b) molecule of symmetry ^E, in which 
case the electronic spin, hence the magnetic dipole mo- 
ment, is completely decoupled from the molecular axis. 
Thus, the electric and magnetic dipoles can indepen- 
dently follow the external electric and magnetic fields and 
the net energy is the sum of the two energies. The be- 
havior of such a molecule in a magneto-electrostatic trap 
does not require special attention and will not exhibit the 
richness of behavior enjoyed by case (a) molecules. 

For real molecules, the dipoles are only approximately 
independent. The spin and rotational degrees of freedom 
are coupled by a Hamiltonian of the form jN ■ S, where 
the electronic orbital angular momentum and the rota- 
tional angular momentum of the nuclei couple to form 
N, with projection A on the internuclear axis (0 in this 
case). couples to the electronic spin S to give total 
angular momentum J. In order to estimate its effect we 
can introduce the spin-rotation interaction as a pertur- 
bation. We have again built a simple analytical model 
including only two rotational states TV = 0, 1. In this 
case, the electric field is the one which lies along the Z 
axis, and we can restrict ourselves to the Mjv = sub- 
set. We assume that 7, J/g, and Um are all much smaller 
than the rotational spacing B. With this perturbation, 
the energy of the ground weak field-seeking state varies 
as oc j^U^Um sin^(Q:)/i?* with the angle a. 



where Ek = Ek{r) are the local energies and |^fe(g')) = 
|0fe(g;r)) the corresponding eigenstates. 

More generally, the internal state of the molecule is a 
superposition of these adiabatic states: 

Mq,t)) = ^afc(t)e--'=W|</.fe(g;r(0)), 



with ujk{t) = J dt' Ek{t')/h. By substituting the ex- 
pansion of |5'((j,t)) into the time-dependent Schrodinger 
equation and solving for the coefficients, we find a system 
of coupled differential equations for the coefficients ak, 



daj 
'dt' 



l-l 



(5) 



From normalization constraints, {4'j\'^\'t'j) = 0, 5jk is 
the Kronecker delta function, and h is Planck's constant. 

We are concerned with transitions from "trapped" 
states \4)t)i which possess potential minima in the cen- 
ter of the trap, to "untrapped" states [(/)„), with maxima 
that push the molecules far away from the trap center. 
The probability pt^u of undergoing a nonadiabatic tran- 
sition during a process can be found to be at most the 
order of [18] 



Vt^u ^ rnax 



\Et-Eu\/h 



and accordingly, adiabaticity requires 



(6) 



(7) 



To avoid calculating time derivatives of the states — 
usually a difficult task — we can recast the matrix ele- 
ments using the Hellmann-Feynman theorem [19, 20], 



h{Et - Eu) ■ 



III. CONDITIONS FOR NONADIABATIC 
TRANSITIONS TO UNTRAPPED STATES 

We turn now to the nonadiabatic changes of internal 

state that lead to trap loss as a molecule moves through 
a trap minimum. Assuming an adiabatic separation in 
the spirit of the Born-Oppenheimer approximation, we 
can distinguish the "slow" center of mass coordinates of 
the molecule, r, from the "fast" internal coordinates (de- 
noted collectively as q). The center-of-mass coordinates 
can be treated classically, since typical translational en- 
ergies are much bigger than the trap energy level spacing. 
Thus, the molecules arc described by position and mo- 
mentum coordinates [r{t),p{t)\. Adiabatic energy states 
are then defined parametrically at each r: 

H{r)\cj>k{q;r)) = Ek{r)\(t>k{q;r)), 



Furthermore, the time derivative can be re-expressed in 
a more intuitive way using the velocity v = ^ oi the 
molecule and = Wr ■ v. The adiabaticity criterion 
then becomes 



\{<i>t\{V rHMu) ■ V\ < 



{Et - Euf 



(8) 



Evaluating Vi?r is relatively simple once the distribution 
of fields is specified. Distinguishing between the different 

radial, polar, and azimuthal spherical components of the 
velocity will be useful, and we can use the expansion: 
V = Vrf + ve9 + v^(j), with f, 0, and (j) the unit vectors in 
radial, polar, and azimuthal directions, respectively. 

As initial examples, we first examine the limiting cases 
of traps with either 11^ = or Um = to gain intuition 
before analyzing traps formed from crossed electric and 
magnetic fields. 
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A. Magnetostatic traps 

In the case of the magnetic quadrupole trap, Eq. (8) 
reduces to the results well-known from previous treat- 
ments [2] . We demonstrate this assuming a magnetic field 
distribution B{r) = VB{x, y, —2z), where VS is the mag- 
netic field gradient on the XY plane. We first consider 
the situation in which particle's trajectory is restricted 
to the XY-plane. Nonadiabatic transitions can con- 
nect the weak field-seeking state \(j)t) = \f,M = -t'f/2) 
to the strong field-seeking state with the same parity 
= I/, M = —1/2). (The first quantum number is 
the parity, and M denotes the projection on the local 
field axis throughout the remainder of this paper.) This 
transition probability is given by 



eral, the dependence on the polar position is given by 



i(*i>.)i = |^. 



(9) 



and the energy gap is \Et — Eu\ ~ Wm- As the magnetic 
field in this plane is given by |K(r)| = VSr, Eq. (7) 
establishes the adiabaticity criterion 



2VBriJ,Bgr, 
2r h 



(10) 



for an azimuthal velocity at a distance r from the 
center. 

One can define a "zone of death" within which the 
molecule will be lost with high probability due to nona- 
diabatic transitions. We can roughly determine its size 
using the radius at which Eq. (10) becomes an equality: 



r0 « {hv^/Ag,npLBVBf'^. 



(11) 



The subscript is a reminder that these transitions occur 
due to the molecules' motion in the <f) direction, orthog- 
onal to the local magnetic field. Spin-flips to untrapped 
states (Majorana transitions) are thus expected within a 
zone close to the center of the quadrupole distribution 
defined by r < r^. In this zone, the change in the ori- 
entation of the field (and thus in the adiabatic basis) is 
rapid when the molecule is on a trajectory close to the 
field zero, and the energy gap between trapped and un- 
trapped states becomes small. Both circumstances being 
essential in Eq. (7), a quick change in the internal state 
may happen. Motion in the radial direction does not 
cause any transitions because the orientation of the field 
does not change along this direction [29]. Thus, the 
size of the zone depends on the azimuthal component of 
the velocity. The estimate given by Eq. (11) can also be 
obtained intuitively, by requiring the Larmor precession 
frequency (proportional to the energy gap) be less than 
the angular velocity (proportional to the rate of change 
of field orientation) , which ensures that the dipole is able 
to adapt to the local field [5, 21]. 

For molecules moving in the XZ- or KZ-planes, the 
radius of the nonadiabatic zone shows the anisotropy in- 
duced by the different gradient along the Z axis. In gen- 



hve 



1/2 



(12) 



.25,„/iBV6(l-h3cos2 6')3/2^ 
Moving across the X- and Z-axes, this reduces to 

4 « {hve/2g^^iBVBY'^, (13) 



and 



re 



{hvg/lQg^^iBVBy^. 



(14) 



The former differs from Eq. (11) due to the component 
of the velocity vq considered. The rate of change of the 
orientation of the field that the molecule experiences nat- 
urally depends on the direction of movement: the change 
is slower when moving in the azimuthal direction (on the 
XY plane, Eq. (11)) than in the polar direction, (on the 
XZ plane, Eq. (14)). In summary, a different radius 
for the zone of death can be assigned to every direction, 
to every nonadiabatic transition from trapped to non- 
trapped state, and to every component of the velocity. 
For the sake of making order-of-magnitude estimates of 
trap loss, however, one may simply use Eq. (11) for every 
position and direction [5] . Finally, it is worth noting that 
Brink et al. [22] have recently calculated loss rates for 
atoms with hyperfine structure trapped in typical mag- 
netic configurations. 



B. Electrostatic traps 

For a quadrupole electric field distribution, £(r) = 
VS{x,y,—2z), with V£ the field gradient on the XY 
plane, nonadiabatic transitions require more elaboration. 
There exist only two different eigenvalues and given the 
degeneracy we have to define a particular adiabatic ba- 
sis. We will use the M projection on the local field to la- 
bel the vectors within each doubly-degenerate subspace, 
adding a subindex 1 or 2 to distinguish M = +1/2 from 
M = —1/2, respectively. Specifically, the molecules are 
trapped in the electric weak-field-seeking states \(pti) and 
1^42) which can make nonadiabatic transitions to the un- 
trapped states \(j)ui) and \(f>u2), which are also degenerate. 
The latter are lower in energy by A from the trapped 
states at zero field (see Fig. 1). 

Using the formalism presented above, we identify two 
types of transitions, spin conserving and spin changing, 
which turn out to depend on different components of ve- 
locity: 



|(<^U|||</'„1)| = |(</'*2|^|</'„2)| = 



(15) 



(A/2)geMe^^rVg 2V(A/2)2 + t/| 

2[(A/2)2 + C/2] ^ n 



mi\j^\M\ = m2\j^\<i>ui^ 



(16) 



ge^leV^V£ ^ 2V(A/2)2 + t/| 



2 V(A/2)2 + C/| 
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Note that they coincide when evahiated at the center of 
the trap, where the symmetry blurs the difference be- 
tween radial and azimuthal velocities. Solving for the 
radius of the zone of nonadiabatic transitions yields for 
the M-changing transitions: 



{hgefXeV<^V£/4: - (A/2)^) 



1/2 



(17) 



This process is completely analogous to the change of M 
in the magnetic field case, modified to accommodate the 
A-doubling. Indeed, if A = 0, then Eq. (17) is identical 
to that for the magnetic case, with magnetic quantities 
replaced by electric ones. Motion in the azimuthal direc- 
tion must be slow enough for the molecule to be able to 
change its dipole orientation with the local electric field. 
Otherwise, it will undergo a nonadiabatic transition. 

Transitions which change only the parity composition 
but preserve M depend on the radial velocity Vr- This 
criterion leads to a slightly different radius for nonadi- 
abatic transitions than that for the M-changing transi- 
tions: 



[(%eMe^.VgA/8)V3-(A/2)^]^/' 



(18) 



This represents a new kind of transition, foreign to the 
usual Majorana loss mechanism involving spins in a mag- 
netic field. For motion in the radial direction, Vr must be 
slow enough for the molecule to change its polarization 
adiabatically with the local electric field. In other words, 
the mixing of parity states / and e that compose the en- 
ergy eigenstate must evolve because a molecule undergo- 
ing radial motion experiences a change in the magnitude 
of the electric field generating the quadrupole trap. 

Notice that the A-doubling plays a vital role in the 
ability of the molecule to undergo a nonadiabatic tran- 
sition. Without lambda doubling, the eigenstates would 
be states of good Q for all field strengths, insuring the 
preservation of these states as the molecule moves ra- 
dially. Moreover, nonzero A-doubling smoothes the po- 
tentials near the avoided crossing, making nonadiabatic 
transitions less likely. Indeed, for < JS? /hgeHeWE 
or, cquivalently, a value of A larger than ^yT^Eg^JIeVS , 
nonadiabatic transitions become very unlikely. A similar 
expression holds for Vr- 

To highlight the difference between nonadiabatic tran- 
sitions in magnetic versus electric quadrupole fields, we 
consider the OH radical as discussed in Ref. [7]. Us- 
ing typical field strengths and molecular velocities of 10 
m/s, and substituting relevant OH (^113/2) energies in 
our J = 1/2 model, the radius of the zone of nonadia- 
baticity would be ~1 fxm for a magnetic quadrupole trap. 
However, such a zone does not exist at all for the electric 
quadrupole trap. Although not necessarily the case for 
all molecules, zones of death for OH in an electrostatic 
trap seem to be essentially negligible at the slow veloc- 
ities obtainable at the terminus of a Stark decelerator. 
The A-doubling effect eliminates this loss channel. 



Loss from the electrostatic trap can also be discussed 
for NH in the metastable state ^A. In this case, the A- 
doublet is significantly smaller than in OH, leading to 
a larger zone of nonadiabatic transitions. We estimate 
that this zone is on the order of 1 /im in diameter for 
an electric field gradient of 3 x 10** V/cm^, which is typ- 
ical for the types of experiments described in Ref. [23]. 
The corresponding loss rate can then be estimated as 
above, yielding a trap lifetime on the order of 10'* s at a 
temperature of 100 mK. This value is much larger than 
the estimated lifetime corresponding to blackbody radi- 
ation [24], and also consistent with the lifetimes found 
experimentally [23] . However, the lifetime decreases to 1 
s when the temperature is 1 mK, and nonadiabatic loss 
processes will become comparable to others upon cooling 
to temperatures of interest. 



IV. TRAPS WITH CROSSED £ AND B FIELDS 

Upon combining £ and 15 fields, both the physics of 
nonadiabatic transitions and its interpretation become 
more complicated. Because the fields are non-parallel 
in general, few quantum numbers are conserved, and 
transitions may occur from the upper trapped state to 
any of the three others in our model. We examine two 
easily-realized experimental configurations. Many more 
are possible, but we start here with the simplest, non- 
trivial examples. Namely, we will consider the influence 
of a uniform electric field on a magnetic quadrupole trap 
as well as the influence of a uniform magnetic fleld on an 
electric quadrupole trap. Nonadiabatic transition rates 
for these situations are found in various limits and loss 
rates estimated. 



A. Magnetic quadrupole Held plus constant electric 

fleld 

We flrst explore the consequences of crossed flelds for a 
magnetic quadrupole trap with a constant electric field, 
as was the situation in the recent JILA experiment de- 
scribed in Ref. [7]. 



1. Trap potential 

In the case of an azimuthally symmetric quadrupole 
magnetic trap, the field distribution near the trap's cen- 
ter is given by S(r) = VB{x,y,—2z). For a simple 
spin-1/2 particle (i.e., ignoring lambda doubling), this 
results in two surfaces, one for weak-field-seeking states 
and one for strong-field-seeking states, which touch only 
in the center where |B| = 0. It is around this region that 
Majorana transitions are expected, as discussed in Sec- 
tion III A. In the case of a molecule with a lambda dou- 
blet, however, there are two such sets of surfaces, offset 
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FIG. 3: (Color online) a) Adiabatic trap potentials for a 
magneto-electrostatic trap in which a constant electric field 
is superimposed on to a magnetic quadrupole field. Slices 
through the center of the trap {z — 0) are shown. They cor- 
respond to the states |1), |2), |3), and |4) which are defined in 
the text. These eigenenergies have been calculated using the 
four-state model (see text) for a molecule with A « 80 mK, 
/ie = 1 D, and /irn ~ 1 /xs (like OH), within a field distribution 
with parameters given by VS — 10'* G/cm, \£o\ — 1.2 x 10'* 
V/cm. Panel b) shows contours of the highest surface |4) 
around the trap center. The arrow points parallel to the elec- 
tric field. 



by A from one another. They can be visualized by rotat- 
ing the four energy curves on the lower panel of Fig. 1 
about its vertical axis. As can be also concluded from the 
figure, the high-field-seeking surface of / parity and the 



low-field-seeking surface of e parity become degenerate 
for a particular value of the magnetic field magnitude. 

Introducing a constant electric field £q = {\So\,0,0) to 
the trap shifts the energies away from the simple, purely 
magnetostatic picture. Figure 3(a) shows the four adia- 
batic surfaces for our model OH molecule for WB =10'' 
G/cm and \£a\ = 1.2 x 10'' V/cm, with eigenenergies 
given by Eq. (4). They are denoted by |1), |2), |3), and 
|4), using the name of the corresponding eigenvectors. 
These surfaces are depicted as slices through the center 
of the trap where z = 0. The upper two surfaces and 
the lower two remain degenerate at {x,y) = (0,0), with 
energies given by the Stark eigenvalues 

±Es = ±y'(A/2)2 + (19) 

where Ueo = ge^J■e\£o\ (a subindex "0" is added to em- 
phasize that the electric field is constant in space). The 
electric field breaks the cylindrical symmetry of the pure 
magnetic quadrupole trap. Thus, the highest energy sur- 
face, |4), possesses roughly elliptical contours of energy 
(Figure 3(b)). Note that the trap is not harmonic, but 
linear in the center. 

The second- highest energy surface, |3), which is also 
able to trap molecules, has a double-well structure due 
to two conical crossings with the next surface, |2), and 
which occur at the coordinates x — ±Es/{gml^B^B). 
At these particular points, electric and magnetic fields 
are collinear and the surfaces can cross. The crossing 
is avoided for any other angular configuration, which 
leads to the conical structure. Thus, while an interesting 
trap geometry may be generated, nonadiabatic transi- 
tions may be possible at these type of points. However, 
we only discuss loss in the central region of the magnetic 
quadrupole trap. 



2. Losses 

As a result of adding the constant electric field to the 
magnetic quadrupole field, nonadiabatic transitions may 
occur between |4) and |3) as in the purely magnetic case, 
but now transitions from |4) — > |2) and |4) |1) are 
possible due to the electric field. 

We begin with |4) ^ |1). Using Ei — -£'4, the condi- 
tion for adiabaticity when moving on the XY-plane (with 
the electric field along the X-axis) is: 



g„ilJ-B^BUeO 



+ Vr sin ( 



< 4 



(20) 



This criterion is always satisfied in the limit Ueo 0, which reduces to the pure magnetic quadrupole field case 
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FIG. 4: (Color online) a) Adiabatic trap potentials for a 
magneto-electrostatic trap in which a constant magnetic field 
is superimposed on to a electric quadrupole field. Slices 
through the center of the trap {z = 0) are shown. They cor- 
respond to the states |2), |3), and |4) which are defined 
in the text. These eigenenergies have been calculated using 
the four states model (see text) for a molecule with A « 80 
mK, — 1 D, and ^rn ~ 1 ^J^B (like OH), within a field 
distribution with parameters given by V£ = 5 X 10=* V/cm^ 
|Bo| = 10'^ G. Panel b) shows contours of the highest surface 
|4) around the trap center. The arrow points parallel to the 
magnetic field 



where there is no such transition. To gain insight from 
this expression, we can assume that the effects from the 
magnetic field are small near the center of the trap ([/„ ~ 
0). This results in the simple expression 



fi9m^J■B^BUeO (W0 cos (j) + Vr siu (/)) ^ 4i?| 



(21) 



The combination (w^ cos (j) + Vr sin (/)) can be easily iden- 
tified as the component (vy) of the velocity that is per- 
pendicular to the electric field, which is along X. Thus, 
Eq. (21) can be re-expressed as 



hgrafJ^B^BVylleO 



< 1, 



(22) 



which may be numerically evaluated in a particular trap 



to check if this transition leads to losses. In fact, the max- 
imum of the left-hand side of the previous expression as 
a function of Ueo is given by (%m^BVi3i'j;)/(3^/^A^), 
which turns out to be neglectable for realistic traps con- 
taining OH. 

Analytical expressions analogous to Eq. (20) can be ob- 
tained for the other two transitions, but their exact form 
is complicated and we do not present them here. In- 
stead, since these transitions are expected to occur near 
the center of the trap, we can assume that the magnetic 
field interaction there is weaker than the electric field 
interaction. We thus consider the magnetic field as a 
perturbation. The degeneracy of the eigenvalues makes 
it convenient to use Brillouin-Wigner perturbation the- 
ory [25, 26] and a zero-order basis which diagonalizes the 
Zeeman perturbation on each parity subset separately. 
For example, Eq. (21) coincides with the expression ob- 
tained using perturbation theory when keeping only lead- 
ing terms for each member of Eq. (8). For the perturba- 
tive series to converge, the condition Um <C Es has to be 
fulfilled. 

Using the perturbation approach to estimate the tran- 
sition rate for |4) |2), which is azimuthal, we find it 
is less likely to occur than the transitions |4) |1) be- 
cause the transition probability is of higher order while 
the energy gap remains nearly the same. Therefore, con- 
dition 22 should set a bound on these losses as well. 

The remaining transition, |4) |3), is the only one 
possible in the absence of an electric field (see Eq. (11)). 
Perturbation theory shows that the previously circular 
zone of nonadiablatic transitions now becomes elliptical. 
The latter is given in polar coordinates by 



r0((/>) 



[(A/2)2 + C/2„] 



Ag^^iBVB [(A/2)2 + U^^ cos^ </,]3/2 



1/2 



(23) 



Note that this reduces to Eq. (11) in the absence of an 
electric field. Substitution of (p — 0,7r/2 in the previ- 
ous expression allows for the calculation of the semi-axes 
of the ellipse (i.e., along the X- and F-axes). The cir- 
cular zone of nonadiabatic transitions shrinks along the 
external electric field and extends in the perpendicular 
direction. 

Assuming that the Z-axis (also being perpendicular to 
the AT-axis) radius of death is analogous to (except 
for a scahng factor due to the doubled gradient), we can 
estimate the scaling of the total losses as a function of 
the electric field magnitude. From the previous expres- 
sions, we see that losses occur within a completely asym- 
metric ellipsoidal zone. The flux of molecules through 
this ellipsoid per time unit constitutes an estimate of 
the loss rate. For low values of the electric field mag- 
nitude UeQ ^ (A/2), the eccentricity of the ellipse in the 
XY plane is given by e = (3/2)i/2(C/eo/(A/2)). If we 
approximate the ellipsoid and the molecular cloud with 
spheres [5], then losses due to the transition |4) |3) 
would be of the same order of magnitude as in the pure 
magnetic quadrupole trap, with corrections on the order 

of (C/eo/(A/2))2. 
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In conclusion, for electric field magnitudes satisfying 
both Ueo (^/2) and Eq. (22), the addition of a con- 
stant electric field to the trap does not increase the rate 
of losses due to nonadiabatic transitions beyond that ex- 
pected from a magnetic quadrupole trap. 

B. Electric quadrupole field plus constant magnetic 

field 

The final crossed-field configuration we will examine 

is that of a electric quadrupole field combined with a 
constant magnetic field. We examine the losses in this 
situation in the same manner described above. 



1. Trap potential 

We again consider an electric quadrupole field con- 
figuration with field distribution £{r) = \/£{x,y,—2z), 
where the gradient V£ is constant. In the absence of a 
magnetic field, there are only two (doubly degenerate) 
distinct surfaces, one weak- field-seeking and one strong- 
field-seeking, and the A-doublet prevents a connection 
between the surfaces in the center of the trap. They 
can be visualized by rotating the two energy curves in 

I 



For a general 4>, which is the angle that the electric field 
makes with the constant magnetic field (assiimed to lie 
along the X-axis), this relation converges to the analo- 
gous expressions for the pure electric trap (see Eq. ( 15) 
and (16)) when Umo goes to zero (wc do not present 
the details of the necessary transformations here, com- 
plicated by a location-dependent change of basis). 

In the case of the azimuthal term in Eq. (24) being 
alone, the ratio of the energy gap to the transition proba- 
bility, thus the adiabaticity, would increase at every point 
with the addition of a magnetic field, but the contribu- 
tion of the radial term seems to depend on the distance. 
The analysis for points very close to the center of the trap 
{Ue -C A/2) shows that the ratio between the left- and 
right-hand sides of Eq. (24) scales as {l + 2Umo/A)~'^. As 
a result, the nonadiabatic transition is actually hindered 
by increasing the magnetic field. 

Expressions analogous to Eq. (24) can be obtained for 
the other two transitions (|4) |3) and |4) |2)), but 
the lengthy analysis will not be presented here. We there- 
fore apply Brillouin-Wigner perturbation theory again to 
consider the electric quadrupole field as a perturbation 
to the homogeneous magnetic field near the trap cen- 



Fig. 1(b). The primary eflFect of adding a small constant 
magnetic field, Bq = {\Bo\,0,0), is to break the degener- 
acy. Higher magnetic fields can modify the connections 
between different potential energy surfaces. For exam- 
ple, tuning the magnetic field can bring regions of sur- 
faces closer to degeneracy, which can induce linear Stark 
dependence along the directions perpendicular to Bq, as 
discussed in Section II A, or induce conical intersections 
similar to those in Fig. 3. 

The surfaces in Fig. 4 correspond to our model OH 
molecule and are generated for fields of V£ = 5 x 10^ 
V/cm^, |Bo| = 10^ G. This magnetic field is small enough 
that no crossings occur between the middle two surfaces, 
and the two trapping surfaces are harmonic. The upper 
surface is characterized by the following spring constants: 
= (3eMeV£)V(A/2),'A:j, = (5eMeV£:)V(A/2 + C/„o), 
and kz = (2geMeVf )2/(A/2 -|- Umo), where UmO = 
5toMb|So| (the subindex "0" is added to indicate a con- 
stant magnetic field). 

2. Losses 

Exact expressions can be found for the probabilities of 
nonadiabatic transitions out of the trapping state, |4). 
The adiabaticity criterion for transitions |4) — > |1) is: 



< 4-^. (24) 

I 

ter. This perturbation theory shows that the transition 
|4) |2) depends on both the radial and azimuthal ve- 
locities. Close to the trap center, both the transition 
probability and the gap between the potential energy sur- 
faces do not change noticeably due to the addition of the 
constant magnetic field. Dependence with the latter ap- 
pears only in second-order terms in the perturbation, Ue, 
far from the center. Thus losses due to this transition are 
not expected to be very different than the ones in the pure 
electric quadrupole trap. However, numerical estimates 
using the exact expressions show that the presence of the 
magnetic field may inhibit the corresponding spin-flips. 

Finally, we consider the transition |4) |3). The 
magnetic field breaks the degeneracy in the upper sub- 
space and nonzero transition probabilities due to both 
radial and azimuthal motion connect both states. Given 
the confining character of state |3) they do not lead to 
trap losses. Regardless, while the transition probabili- 
ties (which connected already the degenerate states in 
the pure electric quadrupole for the analogous choice of 
basis) remains essentially constant, the gap increases lin- 
early with the addition of the magnetic field. Thus these 
fiips are again hindered by the presence of the additional 



/ (A/2)2-hC/2 \ / /^/2)2 \ 

Vd, cos d> UmO H , +VrSm6 \ UmO H , 

/ ^\ V(^/2)' + C^|cos2</. J V(^/2)' + C^|cos2<^ J 
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field. 

In conclusion, the addition of a constant magnetic field 
to an electric quadrupole decreases trap losses beyond 
those found in a pure electric quadrupole trap. In the 
previous analysis we have assumed tiered surfaces like 
that shown in Fig. 4, where Umo < A/2 (magnetic field 
less than ~1.5 kG for our model OH molecule). Degen- 
eracies and crossings between surfaces |2) and |3) that 
occur at higher magnetic fields might deserve special at- 
tention. 



V. REMARKS ON STARK DECELERATION OF 

OH 

The Stark deceleration technique, an experimental 
scheme to obtain samples of cold polar molecules [13], re- 
lies also on the absence of significant nonadiabatic transi- 
tions when switching from one electric field stage to the 
next. Molecules transferred to the wrong states would 
not experience the proper transverse focusing and slow- 
ing and the molecular packet would disintegrate. We can 
apply the formalism above with only minor modifications 
to analyze the Stark deceleration process, and in partic- 
ular, its application to OH in the ^11^/2 state [27]. 

As the molecule moves along the dccclcrator, it experi- 
ences brisk changes in the local field. As a result, the field 
orientation a and strength \£\ may change, giving rise to 
non- vanishing derivatives da/dt and dUg/dt. Eqs. (15) 
and (16) can be modified to include them in our model: 

|(</'n|||.^„i)| = \{Mf^\M\= (25) 

(A/2)geA^e d\e\ 
2[(A/2)2 + dt 

|(<^n|||</'„2)| = m2\j^\ct>ui)\= (26) 

Ue da 
2^{A/2)^ + Ui~dt' 

The changes in the local field are maximum when the 
voltage is switched between pairs of electrodes. We es- 
timate these based on the Stark decelerator and electro- 
static trap of the Lewandowski group at JILA [28] . Dur- 
ing a switching time of ~1 /US, the modulus of the local 
field at the midpoint of a pair of electrodes decreases 
a factor of 30 while rotating by n/2. Assuming that 
the molecule is essentially motionless during this pro- 
cess, this yields maximum values of da/dt = 10^ rad/s 
and d\£\/dt = 10^^ Vcm~^s~^. Introducing in Eq. (6) 
the derivatives given by Eq. (25) and (26), we find nona- 
diabatic probabilities of the order of 10~^ for OH. For 
NH in the ^A metastable state, another Stark deceler- 
ated molecule of interest [23, 28], the contribution due 
to the change of modulus (calculated using Eq. (26)) is 
several orders of magnitude smaller. Thus, nonadiabatic 
transitions are negligible in the decelerator, a conclusion 
in agreement with the findings of [23]. 



VI. CONCLUSIONS 

The presence of both permanent magnetic and electric 
dipoles in case (a) molecules opens the way to rich trap 
dynamics, intrinsically different than that more com- 
monly found in atomic systems. The experimental ex- 
ploration of their behaviour and the potential control of 
their dynamics, which could be exploited through the use 
of arbitrary superpositions of electric and magnetic fields, 
is a question of great importance to current experimental 
research. In particular, we have shown how the simulta- 
neous presence of both fields can lead to new spatial dis- 
tributions within the trap. Also, while homogeneous elec- 
tric fields open new nonadiabatic connections between 
states in a magnetic quadrupole trap, constant magnetic 
fields can inhibit losses in an electric quadrupole trap. 
In summary, for practical design of magneto-electrostatic 
traps, such as those experimentally explored in Ref. [7], it 
should be noted that crossed electric and magnetic fields 
effect can be used to inhibit previously existing nonadi- 
abatic transitions by increasing the energy gaps between 
surfaces or to enhance those transitions by introducing 
new depolarization pathways. The suitable design of the 
field parameters in magneto-electrostatic traps can re- 
duce nonadiabatic transition loss rates below other lim- 
iting loss mechanisms, and thereby open new avenues for 
explorations of ultracold molecular physics without the 
need for building unwieldy loffe-Pritchard traps. 

The comparative analysis of nonadiabatic transitions 
in pure magnetic and electric distributions has also al- 
lowed the identification of a new type of nonadiabatic 
transition in changing electric fields, foreign to the usual 
Majorana loss mechanism. Namely, nonadiabatic transi- 
tion due to the change of magnitude of the electric field 
and not due to its orientation, and this effect is intrin- 
sically linked to the presence of the avoided crossing in- 
duced by the A-doubling. 

We also remark that Hund's case (b) molecules are 
unaffected by the new nonadiabatic loss mechanisms that 
appear in Hund's case (a) molecules, and therefore may 
be more attractive for certain applications. 
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APPENDIX A: STARK AND ZEEMAN 
INTERACTIONS IN CROSSED FIELDS 

When external £- and B-fields do not form a shared 

quantization axis, the scalar products which define Zee- 
man and Stark interaction terms have to be evaluated for 
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general £{£.j.,£y,£z) electric and B(Bx,By,Bz) magnetic 
vectors. Wc outline how matrix elements of the Hamil- 
tonian in Eq. (1) may be calculated in these cases. Both 
electric and magnetic dipole moments are naturally ex- 
pressed in the molecular frame, whose Z-axis coincides 
with the molecular axis. The electric dipole moment lies 
along the molecular axis while the magnetic one is easily 
calculated using electronic operators that are referenced 
to the molecular frame. Fields, on the other hand, are 
best expressed in the laboratory frame. The evaluation of 
the scalar products which define Stark and Zeeman inter- 
action thus requires the expression of dipoles and fields 
in the same reference system, which we can choose as 
the laboratory reference system. Changing the spherical 
components of the dipole moments from the molecular to 
the laboratory frame involves the use of Wigner matrix 
elements (//,; = Y^k^ll^^k)- 

Using these notions, the Zeeman and Stark interaction 
terms can be expressed as: 



where fi and B represent the dipole and the field, re- 
spectively. As the electric dipole moment lies along the 
intcrnuclear axis, its spherical coordinates in the molec- 
ular frame can be expressed as jUe = (0, 0, /Xe). If the 
contributions of the out-of-axis components of the mag- 
netic dipole moment are negligible, as we have assumed in 
the text, Hm « (0, 0, —ijlb{Lo + gSo)/h) is also obtained. 
Thus, the interaction term Al reduces to 



(A2) 



(Al) 



Matrix elements for the interaction operators can be eval- 
uated using the following expression, in the basis without 
parity: 



{JMn\Dl^\J'M'n') = \/2 J + l\/2 J' + 1 J 



J 1 J' 

M q M' 



J 1 J' 



\M-Q 



-n n' 

which, after some algebra, allows one to obtain the following elements in the parity defined basis: 



(A3) 



{JMfl€\D'g*oiie\J'M'fl'e') 



and 



{JM^e\DllnB{Lo + geSo)/h\J'M'n'e') 



-n I n' ) (-i)"-^Ve(A4) 



(A5) 



V2J + W2J' + 1 



J 1 J' 

-M q M' 



J 

-n 



J' 



r 



where the relation = A -|- E allows one to choose the 
correct signs of A and a components in the previous ex- 
pression. The difference in the parity factor in both ex- 
pressions is related to the magnetic dipole moment not 
being a constant (as is fie), but an operator. The evalua- 
tion of the latter — Lo+ggSo — ^in the parity basis extracts 
different signs when acting on the different signed CI com- 
ponents. 



APPENDIX B: SEMICLASSICAL ANALOG FOR 
HUND'S CASE (A) MOLECULES 

The competition between the electric and magnetic ef- 
fects that determines the potential energy surfaces can 
be illustrated by a simple analogy that provides useful 
intuition. The electric dipole moment lies parallel to the 
molecular axis. For a case (a) molecule, the magnetic 
dipole is also constrained to lie along this axis, on av- 
erage, due to strong spin-orbit coupling. Each field will 
thus try to orient the molecule along itself. 

It is reasonable to imagine that the system will essen- 
tially align to the stronger field. That being the case, 
when rotating the smaller field by varying a and hold- 
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ing the magnitudes of both fields fixed, the molecular 
axis will remain anchored to the stronger field while the 
interaction with the smaller field would vary as the pro- 
jection of the dipole on the smaller field, cos a. In fact, 
dependences of the form Um ± Ue cos a and Ue ± Um cos a 
can be extracted from Eq. (4) in the limit of very small 
A. 




FIG. 5: (Color online) "Scissors" model for the polarizabil- 
ity dependence of the electric dipole with the electric field 
magnitude. The molecule works as a pair of identical dipoles 
joined by a ficitious spring that acts to keep them antipar- 
allel. As the field is turned on, both electric dipoles tend to 
align with the field but are partially counterbalanced by the 
spring. The result is a variable electric dipole moment whose 
value depends on the magnitude of the field. 



A semiclassical model can also be found for the case in 
which both fields exert a similar force on the molecule. 
As illustrated in the righthand drawing of Figure 2, if 
the molecular axis makes an angle 71 with respect to £, 
then it makes an angle a — 71 with respect to B, and the 



resulting energy as a function of orientation is 

U{'-fi)±=Ue cos(7i) ± U,n cos(a ~ 71). 

Minimizing (maximizing) this energy with respect to 71 
provides the energies given in Eq. (4) in the absence of 
A-doubling. The ± acounts for the relative orientation 
of the electric and magnetic dipoles. 

To include the effect of A-doubling, we imagine that 
the molecule consists of two, instead of one, antiparallel 
electric dipoles in zero electric field. This corresponds 
to the fact that the molecule in zero field experiences 
no net dipole moment. As the field is turned on, both 
electric dipoles tend to align with the field, but are coun- 
terbalanced by a fictitious spring that acts to keep them 
antiparallel (see Fig. 5). The spring serves to reproduce 
the A-doublet in this model. If the angle between the two 
electric dipoles is denoted 72, then this "scissors" model 
can be described by a classical energy 

f^±(7i: 72) = Ue sin(72/2) cos(7i) ± [/,„ cos(a - 71) 

-A/2cos(72/2), 

(Bl) 

where we have introduced particular dependences on 72 
of both the molecule electric dipole moment and the en- 
ergy of the imaginary spring. The extrema of this func- 
tion constitute a good approximation to the quantum 
energies except that they cannot reproduce the avoided 
crossing. 
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